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Abstract 

We suggest that the 2 + 1 dimensional Gross-Neveu model can give an effective field 
theory description of low-energy quasiparticles in high temperature superconductors. 
The magnetic catalysis of dynamical symmetry breaking is examined. The model shows 
that a magnetic field can induce a phase transition associated with parity breaking. In 
particular, it is intended to give an explanation of a second phase transition observed 
in a recent experiment by Krishana et al. 
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1 Introduction 

The catalysis of chiral symmetry breaking by a magnetic field was studied in Nambu-Jona- 
Lasino models in 2 + 1 jj], 0, ||| and 3 + 10 dimensions. It was also extended to the case 
of QED, external non-Abelian chromomagnetic fields and finite temperatures, curved space, 
super symmetry, etc. || , confirming the universality of the phenomenon. By far it is already 
well known that a magnetic field is a strong catalyst of dynamical symmetry breaking, leading 
to generating a fermion dynamical mass even at the weakest attractive interactions between 
fermions. 

In this paper we promote the hypothesis that this magnetic catalysis of dynamical sym- 
metry breaking can be used to explain the results of a recent experiment on Bi 2 Sr2CaCu 2 8 
(BSCCO) by Krishana et al. ||. At temperature well below the superconducting critical 
temperature (T c ), they discovered a kink-like feature of thermal conductivity at a transition 
magnetic field, which depends on temperature approximately in a quadratic form. They sug- 
gested that the field induced a second phase transition in the superconducting state: a large 
gap appeared abruptly and time reversal symmetry was broken. By now, several possible 
mechanisms of this phenomenon have been discussed |7|, ||, ||, |10j . Laughlin J7| argued that a 



small parity-violating d xy superconducting order parameter was developed, and his hypoth- 
esis led to a model free-energy functional which could display a phase transition induced by 
a magnetic field. But he did not tell how the d xy order parameter was developed. 

In the following, we shall argue that an extended version of the Gross-Neveu model 
|TTL O, ITBf in 2 + 1 dimensions is suitable for studying the low energy excitations of high 



temperature superconductors. The theory shows that an external magnetic field B can 



induce a continuous phase transition at temperatures below the superconducting critical 
temperature: an additional gap is developed and parity is broken. The temperature de- 
pendence of the critical magnetic field shall be calculated and fitted into the experiment 
of Krishana et aZ.f| Unlike those parity-conserving gauge models studied in Refs. || fUj, 
parity is broken in our model. Since parity violation is becoming popular in studying high 
T c superconductors (see e.g. |14|, |7|, [15], [16]]), our model is thus of great interest. 



2 The model 



We limit our discussion to clean two spatial dimensional superconductors in respect of the 
layered structure of cuprates. Following Lee and Wen [17|, we assume that the elementary 
excitations in the superconducting state are well defined fermionic quasiparticles with dis- 



persion E(p) = a/ (e(p) - fi) 2 + A 2 (p), where p = (pi,P2), e(p) = 2t/[cos(pia) + cos(p 2 a)} 
and A(p) = |A [cos(pia) — cos(p2a)] with a the square lattice constant. Obviously, A(p) 
is a d x 2_ y 2-wsNe gap. The Fermi surface consists of four isolated nodes at p = (±^,±^) 
where the gap is vanishing. In the following, we shall consider the case of chemical potential 
/i«0§|. 

In the BCS formulation of the high T c superconductor, the above quasiparticle spectrum 



CTlPj 

4(-p) 



(i) 



can be obtained from the Bogoliubov-Nambu Hamiltonian |]T9 

, | e(p) — zA(p) 

tf = £(c{(p),c i( - P )) 

\^A*(p) -e(p) J 

where c CT (p) and cj.(p) are the quasi-electron annihilation and creation operators with spin 
indices a =|, |. A(p) is related to A(p) by A = exp(i<fi)A where can be interpreted as 



the phase of the complex superconducting order parameter. (By definition, A(p) = |A(p)| 
is real.) Since our aim is to obtain an effective description at low energies and long length- 
scales, it is sufficient to focus on the gapless modes near the d-wave nodes (known as nodons 
||20||), integrating out those electrons far away in the Brillouin zone. Hence, it is legitimate 
to linearize the quasiparticle Hamiltonian ([!]) and perform a rotation in momentum space. 
In the vicinity of each node, say we have 

e(p) -> v F p' 2l 

(2) 

A(p) -> v A p[, 

where p[ = (1/ V2)(-p 1 +p 2 ), p' 2 = (1/ \f2)(-p 1 -p 2 + n / a), v F = 2y/2t f a and v A = A a/V2. 
In the coordinate space, the resulting continuum theory reads 



H = J2j d2x (4( x )> c i( x ) 



^ v F p 2 -/A(p) 1 



c T (x) 

4(x) 



(3) 



y iA*(p) -v F p 2 J 

where the index r runs over four nodal points, pi = —idi with index i — 1,2, and A(p) = 

The presence of an external constant magnetic B field perpendicular to the superconduct- 
ing plane can be incorporated into our effective Hamiltonian ([|) by introducing an external 
(classical) gauge field A and replacing p — > p ± eA with ± depending on whether the 



quasiparticle is electron- or hole-like [pi] , |22| . Here is a subtle point. In order to maintain 
the Hamiltonian both gauge invariant and hermitian, we should not introduce gauge fields 
through the complex superconducting order parameter (or gap) function. That is, we replace 



e(p) — > e(p±eA) whereas keeping A(p) unchanged [E3]. Alternatively, one can understand 



this by an intuitive argument: a gap function A is by definition expressed in terms of the 



paring fields and the interaction potential and therefore should be naturally interpreted as 



a function of coordinate variable x rather than momentum operator p, while e(p) can be 



thought as a kinematic term. Thus, 



H = J2 [d 2 x (c|(x),q(x; 



/ 



vf{P2 + eA 2 ) 



-i£(p) 



\ 



c |(x; 



(4) 



^ iA*(p) -v F (p 2 -eA 2 ) J 
(In our notation, the electron has charge — e.) Accordingly, let us choose a gauge: 

A= (Ai,A 2 ) = (0, Bxi). (5) 

The above Hamiltonian is manifestly invariant under a U(l) gauge transformation A(x): 



A: fx) 



A(x) +<9;A(x), 



e- ieA(x) c CT (x), 



(6) 



A -> e" i2eA(x) A. 



If A(p) had been treated in the same way as e(p), we would have got into trouble of either 
violating the U(l) gauge invariance or having a non-hermitian Hamiltonian 



Following Anderson [ Tl ] , we transform away the (j) dependence of the Hamiltonian (|j) by 

\ 



a unitary transformation H' = U HU 1 with 

/ 

U = 



\ 



e** 
1 



(7) 



/ 



We then have 



H' = / d 2 xxt [< y -iV F (p 2 - eA 2 ) + a 2 v A pi] x r 

r J 

where cr a 's (a = 1, 2, 3) are the usual Pauli matrices and the fermionic fields 



xr=u C\ 



To obtain Eq. fl8|), we have used the fact that 

V0 = -2eA + mv s , (9) 
and assumed that, deep inside a superconductor, the supercurrent velocity v s ~ (i.e., the 



contribution of vortex lattice has been ignored [p4j ). Anderson |22[ suggested that Eq. (g) 
is the effective Hamiltonian appropriate for hole-like excitations in high T c superconductors. 
There exists another unitary matrix U that would lead to an effective Hamiltonian for 



electron-like excitations, identical with (§) except for e — > — e (see Ref. |22| for detail). The 
Hamiltonian (Rp can be equivalently recast into a Dirac ("relativistic") form. To achieve 
this, let us rewrite H' in terms of new fermion fields, related to the spinors Xr via a unitarity 
transformation, 

1 1 /e^ct + ct\ 

* = vf^ +ff3)x ' = viU T -c;); (10) 

and rescale our coordinate space such that xi — > (v&/vd)xi,X2 ~^ (vf/vd)x2, where vd = 
^vpVA takes place of the velocity of light in a true relativistic quantum theory. Also, we 
are free to replace p\ — ► (pi — e^i) for convenience since we are working in the gauge (|5D in 
which A\ = 0. Thus, the effective Hamiltonian for nodal particles (quasiholes, for instance) 
reads 

H ' = VD [ d 2 x^Ml t Di)A, (11) 
where we have introduced the covariant derivative Z), = di — ieAi (with i — 1,2) and the 
2x2 Dirac matrices 7^ = (cri, 02, 03) which satisfy the Clifford algebra {7^,71,} = 25^ u 
{ji,v = 1,2,3). 

We are now in a position to propose a low energy effective field theory for quasiparticles 
in the high T c superconducting state. The partition functional and Lagrangian [F23| read (in 



Euclidean space with h = vd = 1 and (3 = 1/ksT) respectively 

Z = J Vi)V^ e~ ^ dT I d2xC (12) 

and 

C = ^ r ( l3 d T + 7i Di)^ - ^GMv) 2 , (i = 1,2), (13) 

where ij) r = ^73- The 2-component spinor ip r represents the nodal particle (e.g., quasi-hole) 
field with a "flavor" index r = 1,2, • • • , iVf (JVf = 4 for the physical case of four nodes). 
(Hereafter, "flavor" indices r shall be suppressed in our notation.) In this theory, we have 
introduced a four-fermion interaction V\ nt = —^(i J i J ) 2 with g the coupling constant of 
dimensionality of [mass]^ 1 to mimic the quasiparticle interactions in superconducting states. 
The theory ( |T3"D is known to particle physicists as one version of the U(N{) Gross-Neveu 
model in d = 2 + 1 in a constant magnetic field. A mass term is excluded by parity (space- 
reflection invariance), 

ip(x 1 ,x 2 ,T) -> n fiip{-x 1 ,x 2 ,T) . (14) 
As we shall see, it is this discrete symmetry that suffers breakdown. 



Similar models || [T(J were proposed to study the magnetic catalysis phenomenon in high 
T c superconducting states. However, our model is different. To see this, one can rewrite the 
Lagrangian (|13|) as C = ^i{%d T + ^Dj)^ — ^(\I//T\I/;) 2 (where I = 1, 2) by combining two 
2-component fermions in a new four-component spinor and introducing 4x4 matrices 
7 M = diag(cr M , — (Tfj) and r = diag(l, —1). Comparing our theory with those of Refs. || and 



lOfl shows that we have chosen a different four-fermion interaction. Our choice, which shall 
lead to a mass term that breaks parity, is motivated by the popular studies of parity-breaking 
in high T c superconductors [0, 0, [T^ . 



3 Free energy and gap equation 

The quartic term in Eq. flIB|) can be canceled by adding an expression that is quadratic in 
an auxiliary field er(x, r), and that vanishes when a is integrated out. This results in the 
replacement of Eq. (|T3|) with the equivalent Lagrangian 

— — Nt 

^9 

Apparently, C has a stationary point at a = —-^(tpip). Now, the partition functional is 
written as 

Z = j ' V^VipVa e -/<W^f (16) 
To leading order in 1/iVf |2^| , the effective action T[a] is defined by 

e -rH = J v~4)V^e-^ dT $ d2xC . (17) 

The Lagrangian ( p!5|) is quadratic in fermion fields and we can formally integrate out them. 
Thus, 

Nt 

I>] = /?V 2 -V - ln(det K), (18) 
where V% — J d 2 x is the area of the space and the "matrix" 

^rxr.r'x'r' = -['jiD i +'j 3 d T + a]5 2 (x-x')5(T-T')5 rr >. (19) 

We only have even powers of a because the trace of an odd number of Dirac matrices 
vanishes (alternatively, because a changes sign under parity transformation (13)). At first 



glance, it would seem that (a) = is automatically a stationary point since T is even in a. 
However, if we are interested in the breakdown of parity, the issue is precisely whether there 



are stationary points other than (a) = 0. Fortunately, to settle this issue, we only need to 
compute T for constant a p?| . 

Since K is invariant under imaginary time translations, it is convenient to perform a 
Fourier transformation on r: 

I*! 3 fin- f 13 rl-r' 

K , , , = / — e iuj " T / —e-^x 
Jo VP JO VP 

= -[jiDi - iuJnl-i + <t]5 2 (x - x')5 nn >5 rr >, (20) 
where uo n = (2n + l)vr//3 is the usual Matsubara frequency for fermions. Then, we have 

/+oo 
d 2 x ^2 tr {(x| ln(7,A - ^n73 + o")|x)} 

n=— oo 

n +0O 

= -N t / d 2 x tr{(x|ln((z 7i A) 2 + ^ + a 2 ) |x)} (21) 

n=— oo 

in virtue of the fact that 

22 tr ln {li A - ^nT3 + tr) = y^trln (73 [7jDj - ^73 + cr] 73) 

n n 

= tr In (-7^ A + ia;_ (n+1) 73 + a) = tr In (- [7; A - K*Vys] + cr ) • (22) 

n n 

Therefore, ln(det K) can be expressed in terms of an integral over the proper time s, 

ln(detK) = --J2 f d 2 x [°° ^tr(x|e- ls [^ A ) 2+ ^ +CT2 ]|x), (23) 
2 „ J Jo s 



where (ryjA) 2 = —DiDi + ieB^i^- Following Schwinger's proper time approach [p8l, we 
find the matrix element 

( x '| e -^(n^i) 2 | x /') = ZlR^y-^ s )eB[cot(eBs) + z 73 ], (24) 

47T 



where 



R(x',x"; s) = exp[-(x' — x") 2 eB cot(eBs) + ie dxiAi(x)], 
4 Jx" 



in which the integration path is a straight line. Substituting Eq. (|24| ) into Eq. (p3|), we get 



ln(det K) 



iN f eB 

47T 



cot(eBs). 



(25) 



Inserting Eq. (^) in Eq. ( p~8|) gives us the effective action. Let T[a] = pVzFfo) where F(a) 
is called the free energy. After a Wick rotation s —is, 



„. , N f a 2 N f eB ^ f°° ds _ 



2# 4vr/5 



-s(ul+a 2 ) 



coth(e-Bs) 



(26) 



where we have explicitly introduced the ultraviolet cutoff (Debye frequency) A. An integral 
very analogous to that of Eq. fl26|) was evaluated in Appendix B of Ref . |l] . We thus skip all 
intermediate analyses and directly go the final result of the free energy, 



Ztt 



V2. 



2?2 



a 



N f 



2vr/?/ 2 



ln(l + e~^) + 2 ln (! + e~^° 2+2k ' 12 ) 



k=l 



+ o( x ) 



(27) 



4= is of the dimension of mass, 

g v% ' 



where I = \eB\~ x l 2 is the magnetic length, Mq 
and C{z,q) is the generalized Riemann Zeta function f29|. This result is exact for arbitrary 
strength of B. 



The gap equation, OF /da = 0, reads 



a 



Ma 



1 ,A a 2 l 2 1N 1 , /3<r 
v^Z v 2 2 ; 2 V 2 



(o-2 + 2Jb/Z 2 )-V2 

+ 72 



:28i 



^ (e/3v ^7P + 1) _ 

The nontrivial solution to Eq. fl28|) yields the quantum thermal average (a). As Eq. (|TJ|) 
implies, a non- vanishing (a) means that fermions gain dynamical mass, m = (a). In other 
words, a gap opens at the nodes of the quasiparticle spectrum. Meanwhile, parity is broken 
by the mass term (a) (ipip) to appear in the effective Lagrangian. As B — > 0, Eq. (p8|) can 



10 



be explicitly solved, yielding the known result |L2 



(a) = -M + 2/T 1 ln(i[l + y/1 - 4exp(/?M )]). (29) 

Therefore, for some range of parameter space (T, Mq) there is spontaneous parity breakdown. 
Obviously, as T — ► 0, Eq. (p9|) has a nontrivial solution only if M < (i.e., <? > 7r 3 / 2 /A). 
For finite -B, the magnetic field however changes the situation dramatically. Firstly, an 
external magnetic field itself explicitly breaks the parity invariance of the Lagrangian ( (13f) 
since it requires that the Lagrangian be invariant under transformation (14) together with 
transformations {A\, A 2 ) — > {—A\, A 2 ) and B — > —B. Secondly, as we learned from Gusynin 
et al. the magnetic field can catalyze parity breaking. This shall be discussed further in 
the following section. 

4 Comparison with experiment 



Next, we shall investigate the free energy fl27D and the gap equation (Pq ) numerically for the 



experiment of Krishana et a/.0. The theory has two free parameters, g and vp. The latter 
determines the typical energy scale of the theory. (In this paragraph only, we shall explicitly 
restore h and vd in order to discuss experiments.) The lattice constant for the cuprate 
BSCCO is a ~ 5.41A(see, for example, |H|]). The point-contact tunneling experiment on 



BSCCO |32| shows that the ratio of 2A //cb^c is roughly 11.6 ~ 12.4. (T c denotes the 



superconducting critical temperature, and A is the gap associated with the occurrence of 
superconductivity and needs not to be confused with the gap given by Eq. Q2"8|).) If we take 
12.4 as the ratio and T c = 92K f|, then A « 49.15meV. This gives hv 2 = A a/V2 w 

11 



0.188eV-A. On the other hand, the superconducting coherence length for BSCCO can be 
taken typically as £o = 20A[|nj. By definition, hv F = 7rA £ ps 3.09eV-A. Therefore, we find 
hv D = U^vpV2 ~ 0.762eV-A. The Debye cutoff A is of the order of Hvd/cl ~ 0.14eV, and 
should be regarded as the highest energy scale in our effective field theory. 

The free energy ([27]) is plotted in Fig. [I]. It exhibits that at low temperatures a magnetic 
field induces a (continuous) phase transition: F(a) can have a stationary point (cr) ^ 0, which 
infers that parity is spontaneously broken. As shown in the insert of Fig. [I], a slight increase 
of the field above the critical point (5B = 0.2T) results in a gap 2(a) ~ 16K at T = 10K. It 
follows that the quasiparticle density is substantially suppressed by a factor exp(— 2/3 (a)) ~ 
20%. Fig. |2] shows the critical line of B-T determined by Eq. ( pB"D in comparison with the 
experimental data from Krishana et al. 0. We find that the experimental data can be well 
fitted by a very small Mq []33|| . According to our theory, the critical line means that parity 
is restored above some temperature that depends on B. We thus conclude that the basic 
feature of the experiment || can be well explained by the present theory. 
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Figure 1: The free energy F' = F(<r) - F(0) as a function of a at T = 10K, M = 0, and 
Hvd = 0.762eV-A. Insert: (a) B = 1.44T is just below the critical point; and (b) B = 1.64T 
slightly above the critical point. 
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Figure 2: The critical line of B—T with parameters M = and Hvp = 0.762eV-A. The 
experimental data is the courtesy of Krishana et al. 
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